Four Wave Mixing from Fe 3+ Spins in Sapphire 
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ions in sapphire exhibit an Electron Spin Resonance (ESR) which interacts strongly with 
high-Q Whispering Gallery (WG) modes at microwave frequencies. We report the first observation 
of a third-order paramagnetic nonlinear susceptibility in such a resonator at cryogenic tempera- 
tures, and the first demonstration of four wave mixing (FWM) using this parametric nonlinearity. 
This observation of an all-microwave nonlinearity is an enabling step towards a host of quantum 
measurement and control applications which utilize spins in solids. 



Since the development of the first laser [T] , a multitude 
of nonlinear effects have been observed in optical systems. 
Optical second- [2] and third-harmonic generation |3J S] , 
optical sum- frequency generation [5], optical parametric 
oscillation and amplification [6] [7], Raman lasing [8], and 
two-photon absorption [9] are all well-characterised non- 
linear effects which have been instrumental in the devel- 
opment of the past few decades of modern optics. High 
quality optical cavities allow the effect of the nonlinear- 
ity to be greatly enhanced, and have lead to many new 
applications including the implementations of frequency 
combs through parametric frequency conversion effects 
[TUJ-[T2]. Optical nonlinearities are crucial for switching 
and modulation in modern communications technology, 
and are an enabling capability for future implementa- 
tions of optical computer technologies, including the pos- 
sibility of a quantum computer based on encoded sin- 
gle photons [15] . Recently, dramatic progress has been 
made in using microwave systems for quantum informa- 
tion and measurement, with nonlinearities playing a crit- 
ical role. Josephson junctions in particular, which op- 
erate at microwave frequencies, act as a nonlinear in- 
ductor which permits uneven spacing of energy levels, 
leading to individual addressability of energy states us- 
ing an external field. This, and other strongly nonlinear 
systems are currently of considerable interest for a new 
generation of quantum measurement experiments includ- 
ing quantum- limited amplification |14j . single quadra- 
ture squeezing with tunable nonlinear Josephson meta- 
materials [T3], readout of superconducting flux qubits 
|16j . and frequency conversion with quantum-limited ef- 
ficiency |17j . An addressable quantum memory with co- 
herence times long enough for quantum computing appli- 
cations could potentially be achieved through the manip- 
ulation of electron spins in a crystal lattice host, which 
typically occurs at microwave frequencies, and can have 
characteristic relaxation times of order seconds. This, 
along with the potential for large collective couplings, 
have provoked great interest in electron spins in solids as 
potential quantum memories for superconducting qubits. 
In particular, nitrogen-vacancy (NV) centers in diamond 



[18], Cr 3+ spins in sapphire [19], and nitrogen spins in 
fullerene cages & phosphorous donors in silicon [3U] have 
been well studied in circuit QED experiments coupling 
superconducting resonators to electron spin ensembles. 

In this Letter, we demonstrate the resonant enhance- 
ment of the weak nonlinear \^ paramagnetic suscepti- 
bility present in a parts-per-billion concentration of elec- 
tron spins in sapphire. To the authors' knowledge, this is 
the first observation of a paramagnetic nonlinear process 
purely at microwave frequencies in a crystalline host. De- 
generate FWM is achieved with the application of only 
a single pump field, with the pump and idler frequencies 
enhanced by ultra-high Q-factor WG mode resonances. 
FWM is an enabling process for both frequency comb 
generation and many quantum computing and metrol- 
ogy applications. Our system is further suited to these 
applications due to the extremely low dielectric loss tan- 
gent at millikelvin temperature, which persists even at 
single photon input power |21j . 

The experimental system consists of a cryogenic sap- 
phire resonator-oscillator p2H5Tj as shown in Fig. [lja). 
The system is cooled to liquid helium temperature and 
pumped with microwave power to excite WG mode res- 
onances. As a result of the manufacturing process, para- 
magnetic Fe 3+ are included in the sapphire lattice at 
a concentration of 150 ppb (~ 10 16 spins in the lat- 
tice) [37] ■ The crystal field splitting results in an inho- 
mogeneously broadened electron spin resonance (ESR) 
with 27 MHz lincwidth [28] at zero applied DC magnetic 
field, corresponding to the spin- 1 1/2), |3/2), and 1 5/2) 
states of the ion. Within the system a complex interac- 
tion occurs between the microwave input field, a dilute 
paramagnetic Fe 3+ spin system, and 27 A1 lattice ions, 
which ultimately results in the production of signal and 
idler photons equally spaced in frequency characteristic 
of degenerate FWM. The resonator-oscillator geometry 
is such that two microwave resonances exist within the 
ESR bandwidth, as shown in Fig. [ijb), which act to res- 
onantly enchance both the pump and idler fields. The 
pump resonance frequency of uiq = 12.0375 GHz is co- 
incident with the 1 1/2) —> 1 3/2) transition residing at the 
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FIG. 1. (a) Cross section of the experimental package. The sapphire resonator is shown in its cavity, which is mounted in 
a vacuum chamber at the end of an insert in a liquid helium dewar. Herotek DT8016 power detectors are used at the output 
to generate a voltage proportional to incident microwave power, measured with an oscilloscope, (b) (Not to scale.) Schematic 
of the system described by our theoretical model. Qo and Q- represent fixed microwave WG mode resonances in the sapphire 
resonator with bandwidths of order 10 Hz. fi + models a lossy resonance at Q + = 2Qq — f2_ as no WG mode exists at 12.045 
GHz. The applied pump frequency ojq can be selected in a range of over 4 kHz around Qo to successfully result in the generation 
of lo- whose frequency only changes over a narrow range <40 Hz, and ui+ with a frequency range of order 8 kHz. 



(iii) 



(iii) 



E 



0) 



(ii) 



Spectrum Analyzer Trao 



(ii) 







Time (sec) 



(ii) 



: (i) 



6 8 10 

Time (sec) 



FIG. 2. Power detected in transmission through the resonator when (a) the 12.037 GHz pump signal is offset 2.671 kHz above 
resonance, and (b) when the 12.037 GHz pump is offset 2.659 kHz above resonance. In both cases, the excitation signal was 
switched on at zero seconds. 



maximum of the ESR, and the idler resonance frequency 
of w_ = 12.0298 GHz is in the wings of the ESR. No WG 
mode exists within the ESR bandwidth at the signal fre- 
quency UJ + . 

At low pump powers, excitation was only observed at 
the pump frequency. However, as is characteristic of de- 
generate four-wave mixing, after a threshold power level 
is surpassed, continuous excitation of the signal and idler 
fields is also found to be present, with the signal and idler 
frequencies equally spaced Aw = 7.669 MHz from the 
pump frequency, and the idler frequency clamped to its 
whispering gallery mode resonance. Remarkably, signal 
and idler excitation was observed to appear anywhere 
from instantaneously to ^7.5 seconds after application 
of the pump, with the delay being strongly dependent on 
the pump power and its detuning from resonance. Fig- 
ure [2] shows the transmitted power through the resonator 
as a function of time for two different cases, with insets 
showing the spectrum analyser trace before (Fig. [2ja)(i)) 
and after applying the synthesizer signal. In the first 
case, when the synthesizer is switched on (Fig. [2Fa)(ii)) 



the transmitted power remains constant for 7.58 seconds, 
after which time the signal and idler appeared simultane- 
ously (Fig. f2fa)(iii)) and the total output power detected 
was seen to rise. In contrast, in the second example the 
signal and idler appear significantly faster, but relax over 
a period of several seconds before reaching a steady-state 
transmitted power level. 

The complex time dynamics of the interaction can be 
qualitatively understood as being due to the slow seeding 
of the parametric process due to energy transfer in the 
Fe 3+ spin system. Upon application of the pump field, 
the sub-set of Fe 3+ ions within the ESR at the pump 
frequency begin to absorb energy. The hyperfine lat- 
tice interaction between the individual Fe 3+ spin packets, 
and the Al nuclear spins then slowly transfers power 
through a cross-relaxation process from the pump mode 
frequency, down to the idler resonance frequency thus 
seeding the signal. A similar behaviour has been ob- 
served in optical systems, where Raman scattering from 
the pump into the signal seeds the parametric process, 
and significantly reduces the threshold power for obser- 
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vation four- wave mixing [251 130] . In the optical case, the 
Raman scattering occurs virtually instantaneously. Here, 
by contrast, the time constant of the hyperfine lattice in- 
teraction can be extremely long, with previous studies 
with Cr 3+ doped sapphire [3T] recording transient relax- 
ation times on the order of 5 seconds. We attribute the 
complex dynamics observed over long time scales in our 
experiment to this fact. 

The four-wave mixing operates for pump frequencies 
over a range spanning 4 kHz, corresponding to ^400 
times the linewidth of the pump WG resonance. Fig- 
ure [3] shows the output power of the signal and idler 
frequencies for a selection of input powers swept over the 
pump WG mode resonance frequency. Four-wave mixing 
is a phase sensitive process, and we observe an apparent 
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FIG. 3. Output power of the signal and idler as the input 
pump is swept in frequency over the WGHi7,o,o resonance. 
The power of the signal is upshifted by the ratio of the am- 
plitudes a_/a+. This upshifted curve nearly directly overlaps 
the signal curve. The transmission curve of the WGHit^o.o 
resonance is shown for reference. 



strong phase-mismatch when the pump frequency tunes 
closely to pump WG mode resonance. A large enough 
phase mismatch ensures that four-wave mixing is effec- 
tively suppressed. Due to the resonant enhancement of 
the idler field, it's frequency is clamped strongly to the 
signal resonance frequency, varying by less than 36 Hz 
over the full pump frequency range as shown in Fig. [4] 
This allows the signal frequency to be widely and pre- 
dictably tuned by tuning the pump frequency. The tun- 
able bandwidth decreases with pump power and is, for 
example, only several hundred Hertz at a pump power of 
5 dBm. 

The full model of the parametric process including 
cross-relaxation induced seeding is beyond the scope of 
this article. Here we instead neglect the seeding process, 
and estimate the parametric nonlinearity through a sim- 
ple three mode picture with Hamiltonian: 
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where the terms on the first line are the rest energy of 
the system, while those on the second account for the 
nonlinear interaction, with g being the nonlinear interac- 
tion strength. The annihilation operator hj describes the 
field amplitude in mode j, with (a\aj) being the mean 
photon number in the mode and the subscripts 0, +, 
and — respectively denoting the pump, signal and idler 
modes. Qj is the resonance frequency of mode j. The 
splitting of the WGHi7,o,o mode at f2o is neglected since 
the splitting frequency is far smaller than the spacing of 
the pump, signal, and idler frequencies, and thus it is 
expected not to contribute significantly to the physics. 

Applying the quantum Langevin Equation to Eqn. [I] 
[32], one can then find equations of motion for the pump, 
idler and signal. In the rotating frame, this yields the 
expectation value equations 
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where aj = (dj), jj and Aj are, respectively, the de- 
cay rate and detuning from resonance of field j, 7o,i n 
and ao.in are the input coupling and amplitude of the 
incident pump field, and since the idler and signal are 
not pumped, a-^n = <%+,in = 0. Expressed in terms of 
half-bandwidths, 7_ = 6 Hz, 70,; = 5 Hz, and 70. u — 
6.7 Hz. Consistent with our experiments, we model the 
signal resonance as a lossy resonance such that the sig- 
nal dynamics are fast compared with the pump and idler 
dynamics. Hence, 7+ represents a lossy damping and is 

a 2 

related to the ratio of the amplitudes 7+ = 7_ . Equa- 

@ a+ . 
can then be adiabatically eliminated, giving an 
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equation of motion for the signal: 
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The effective nonlinearity g' can be related to the intrin- 
sic nonlinearity g, and expressed in terms of only the 
pump parameters, given by: 
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Here, the steady-state intracavity pump amplitude ao is 
related to the threshold power by i* hrcsh = hn \a^ sh \ 2 . 
Figure [5] shows the effective nonlinearity and threshold 
power as a function of normalised detuning from reso- 
nance, calculated using the parameters of our system. 

In summary, we have demonstrated that a strong 
nonlinearity at microwave frequencies, arising from only 
a parts-per-billion concentration of paramagnetic ions, 
leads to degenerate four-wave mixing in a cryogenic 
sapphire resonator when pumped with only a single fre- 
quency. Long characteristic times on the order of several 
seconds were observed due to slow cross-relaxation and 
interaction with lattice ion nuclear spins, and broad tun- 
ability can be achieved by altering the pump frequency. 
Our system has the potential for application in a host of 
future quantum computing and metrology experiments 
where low microwave loss and strong nonlinearity is 
desirable, such as measurement of qubits in circuit QED 
setups, single quadrature amplification and squeezing, 
quantum-limited parametric amplification, or potential 
use in nanoscale magnetometry [33] with the benefit 
of having the necessary amplification and nonlinearity 
integrated within a single device. 
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